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ABSTRACT 

We present the 2 — 8 keV number counts from the 1 Ms Chandra observation of the Chandra Deep 
Field North (CDF-N). We combine these with the number counts from a 78 ks exposure of the Hawaii 
Survey Field SSA22 and with the number counts obtained in independent analyses of the CDF-S and 
the Hawaii Survey Field SSA13 to determine the number counts from 2 x 10^^^ to 10^^^ erg cm^^ s^^. 
Over this flux range the contribution to the X-ray background is 1.1 x 10~^^ erg cm^^ s^^ deg^^. When 
the contributions above 10~^^ erg cm~^ s^^ from BeppoSAX or ASCA observations are included, the 
total rises to 1.4 x 10~^^ erg cm~^ s~^ deg~^. However, there appears to be substantial field-to-field 
variation in the counts in excess of the statistical uncertainties. When the statistical and flux calibration 
uncertainties (in both the background and source measurements) are taken into account, as much as 
0.5 X 10"^"'^ erg cm~^ s~^ deg~^ could still be present in an unresolved component. 

Subject headings: cosmology: observations — galaxies: distances and redshifts — galaxies: evolution — 
galaxies: formation — galaxies: active — galaxies: starburst 



1. INTRODUCTION 

The X-ray background (XRB) was the first cosmic back- 
ground detected (Giacconi et al. 1962) and has been exten- 
sively characterized (Fabian & Barcons 1992). Its photon 
intensity P{E), where E is the photon energy in keV and 
P{E) has units [photons cm~^ s~^ keV""'^ sr"'^], can be 
approximated in the 1 — 15 keV range by P{E) — AE^^ 
where F ~ 1.4 (e.g., Marshall et al. 1980; Gendreau et 
al. 1995; Chen, Fabian, & Gendreau 1997). However, the 
XRB normalization is still somewhat uncertain with values 
in the 2 - 8 keV band lying between 1.3 x 10"" (MarshaU 
et al. 1980) and 1.8 x 10"" erg cm'^ s'^ deg'^ (Vecchi 
et al. 1999). 

After the discovery of the XRB there was controversy 
over whether it arose from a superposition of discrete 
sources or from thermal bremsstrahlung from a hot intcr- 
galactic gas (e.g., Field 1972). We now know that the XRB 
cannot originate in an uniform hot intergalactic medium 
since the absence of a strong Compton distortion in the 
cosmic microwave background spectrum puts a stringent 
upper limit (~ lO"'') on such a contribution (Wright et 
al. 1994). However, there may be other sources of hard 
X-ray photons (e.g., Abazajian, Fuller, & Tucker 2001) 
which could contribute to the XRB without producing a 
Compton distortion, so constraining any residual diffuse 
background is of considerable interest. 

The sources contributing the bulk of the 2 — 8 keV XRB 
were not found prior to the launch of the Chandra X-ray 
Observatory. While the soft (0.5 — 2 keV) background 
was largely resolved into sources by the i? 05*^ T" satellite 
(Hasinger et al. 1998), the spectra were too steep to ac- 
count for the flat XRB spectrum. With the launch of 



Chandra and the XMM-Newton Observatory the situation 
has evolved rapidly. Early 100 — 300 ks Chandra obser- 
vations determined the number counts in the 2 — 8 keV 
band above 10~^^ erg cm~^ s~^ (Mushotzky et al. 2000; 
Giacconi et al. 2001; Tozzi et al. 2001) and resolved the 
majority of the 2 - 8 keV XRB. Hasinger et al. (2001) 
have reported parallel results with XMM-Newton. Now 
the two 1 Ms Chandra exposures of the CDF-N (Brandt 
et al. 2001, hereafter BOl) and the Chandra Deep Field 
South (CDF-S; Giacconi et al. 2002, hereafter G02; Cam- 
pana et al. 2001, hereafter COl) have extended the counts 
to 2 X 10^^^ erg cm^^ s^^. However, there is considerable 
spread (about 40%) in the normalizations, with the counts 
in the SSA13 (Mushotzky et al. 2000) and CDF-N (BOl) 
flelds being higher than the CDF-S counts (G02; COl) and 
the XMM-Newton Lockman Hole counts (Hasinger et al. 
2001). 

Here we provide a more detailed analysis of the number 
counts in the CDF-N, modelling the effects of Eddington 
bias and incompleteness. We also analyze a new 78 ks 
observation of the SSA22 field. We combine our counts 
with those measured by Mushotzky et al. (2000) from the 
100 ks observation of the SSA13 field and those measured 
by COl from the 1 Ms observation of the CDF-S to aver- 
age over a spread of fields. We are currently analyzing the 
latter two fields using the same methods described here 
(L.L. Cowie et al., in preparation), but we note that simi- 
lar number counts have been obtained for the CDF-S using 
very diflFerent methods (G02; COl). 

We break with the tradition of X-ray number counts by 
using differential and not cumulative counts. Differential 
counts have many advantages in the statistical indepen- 
dence of the data points and in the ease with which breaks 
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and shape changes may be seen (c.f., Jauncey 1975). 
2. X-RAY DATA 

The observations of the CDF-N are presented in BOl, 
where the 2 — 8 kcV image and the exposure map are de- 
scribed. We follow BOl and use 2 — 8 keV fluxes since there 
is little sensitivity above 8 keV in the data. When we use 
the data of other authors, we convert to this band using 
consistent assumptions for the spectral index. 

An important issue is the choice of spectrum for con- 
verting counts to flux. BOl converted using the F value 
determined from the ratio of the hard to soft band counts 
for each individual source. Others have converted using a 
single r chosen to match the expected average spectrum 
of the sources. G02 and COl used F = 1.4, corresponding 
to the shape of the XRB, while Mushotzky et al. (2000) 
and Barger et al. (2001) used F = 1.2, based on the ar- 
gument that the average spectra at these fluxes must be 
slightly harder than that of the XRB once allowance is 
made for the contribution of brighter sources (which, on 
average, are softer than the XRB). None of these methods 
is fully satisfactory for individual sources whose spectra 
may be more complex than a simple power-law. We adopt 
the F = 1.2 method which should give the best average 
conversion and hence the best estimate of the contribu- 
tion of the sources to the XRB. This F, together with the 
Nh = 1.6 X 10^° cm~^ of Galactic absorption, gives a con- 
version of counts to flux of 2.19 x 10^^^ erg cm^^ for the 
charge transfer inefficiency corrected images of BOl. The 
F = 1.4 conversion would reduce the fluxes by 4%, while 
use of individual Fs would increase them by 13%. In com- 
bining the new counts with published number counts, we 
have applied the appropriate conversion. 

There are various techniques for source selection, includ- 
ing sliding cells, wavelet based algorithms (BOl; COl), and 
the optimal matched filter method (Vikhlinin et al. 1995). 
In the present analysis we use a matched aperture, which 
consists of a circular aperture of specified diameter and a 
background subtraction region which is an annulus 2 — 3 
times the aperture radius. We determined the background 
from a cleaned image in which all known sources were 
masked out with artificial noise appropriate to the local 
exposure time. The background-subtracted source counts 
were determined by centering apertures on each pixel in 
the image and finding the local maxima. In each case if 
there was a secondary maximum within 2" of the primary, 
we combined the two sources. We determined the average 
background per aperture and selected only sources whose 
background-subtracted counts lay above the level where 
Poisson noise in the average background within the aper- 
tures would provide a probability of one false detection in 
the area. 

We first analyzed the central area of the CDF-N im- 
age lying within 2.5' of the optical axis. Here we used 
a 2" aperture selected to optimize the recovery (based 
on the Monte Carlo simulations described below). This 
gave an average number of background counts of 1.8 per 
aperture; thus, we selected sources with background- 
subtracted counts lying above 9.0. We next analyzed 
the 6.5' radius area around the optical axis (including 
the previously analyzed 2.5' area) using a 4" diameter 
aperture and selected sources with background-subtracted 
counts above 14.7. The regions have fairly uniform ex- 



posure times, which justifies the use of the average back- 
ground. However, to test this, and also to check that 
the interchip regions, where the properties of the sources 
might be unusual, were not affecting the counts, we reran 
the search, clipping out the regions with lower exposure 
times. This had no significant effect on the number counts. 
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Fig. 1. — The ratio of the measured number counts to the input 
number counts in the 2.5' radius (solid boxes) and 6.5' radius (open 
boxes) regions. The width of each bin is shown by the short horizon- 
tal lines. Vertical error bars are smaller than the filled symbols. At 
high fluxes the retrieved counts precisely match the input counts. 
At intermediate fluxes Eddington bias becomes important and the 
retrieved number of sources exceeds the input counts. In the lowest 
flux bins, source recovery becomes incomplete. The corrected area 
is largest in the 6.5' region analysis except in the lowest bin where 
the star shows the recovered fraction in the 6.5' region times the 
relative area of the two regions. 

We used the source counts determined in the 4" diame- 
ter aperture, with the subtracted background determined 
in the 8" — 12" annulus, to get the fluxes of our sam- 
ple. These counts were converted to count rates using the 
BOl exposure map. We determined the enclosed energy 
correction by measuring the ratio of the counts in a 20" 
diameter aperture relative to the 4" diameter aperture for 
all sources with fluxes above 10~^^ erg cm~^ s"'^ and flt- 
ting this with a third order polynomial as a function of the 
off-axis radius. There are no extended sources included in 
this computation. The enclosed energy correction agrees 
well with the standard ACIS PSF correction. 

The flux calibration was tested by comparing the de- 
tected sources with the BOl catalog converted to F = 1.2. 
Above a flux limit of 3.5 x 10~^^ erg cm~^ s~^, the sample 
within the 6.5' radius region contains 152 sources and BOl 
141, of which 134 are common. For the sources in the BOl 
catalog above 5 x 10^"'^*' erg cm^^ s^^, the mean fluxes 
agree to 1% and the dispersion of the fluxes, when multi- 
ple sources are excluded, is 8%. Number counts generated 
from the two catalogs are nearly indistinguishable above 
this flux. At fainter fluxes the catalogs become more dis- 
joint, reflecting the different source selection procedures, 
but the fluxes of the overlapping sample continue to agree. 

The effects of incompleteness and bias were measured 
using Monte Carlo simulations in which a number of model 
X-ray images were created with noise properties matching 
that of the 1 Ms fleld and with sources drawn from pa- 
rameterized descriptions of the number counts. 100 such 
images were generated and analyzed in the same way as the 
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actual X-ray image of the 1 Ms field. The artificial images 
were generated by masking all known sources from the real 
1 Ms image and substituting artificial noise. Sources were 
then randomly drawn from the number counts model and 
added to this artificial image. To match the point spread 
function (PSF) the template sources were drawn from the 
brightest source in the real 1 Ms image in a 1' width an- 
nulus centered on the input position. This template was 
background-subtracted and flux-scaled to produce the cor- 
rect input flux. The source identification and cataloging 
procedures were now run on the artificial image to gener- 
ate the output counts. The output counts averaged over 
the 100 simulations were then compared with the input 
counts to determine the recovered fraction, which could 
then be used to correct the measured number counts. The 
input number counts model was next adjusted to match 
the observed number counts and the process repeated. The 
recovered fractions are shown in Figure |l| and illustrate the 
effect of Eddington bias (the recovered sources are bright- 
ened, which causes the observed counts to be high at in- 
termediate fluxes) and the effect of incompleteness in the 
recovery at the lowest fluxes. 

For fluxes between 2 x 10~^^ and 3x 10~^^ erg cm~^ s~^, 
we used the 2.5' radius sample and its corrected area. 
Between 3 x 10^^^ and 5 x 10^^^ erg cm^^ s^^, we 
used the 6.5' sample and its corrected area. Above 
5 x 10~^^ erg cm~^ s~^, we used the entire CDF-N area 
lying within 10' of the optical axis, an area of 294 square 
arcminutes. At these large off-axis angles, simulations are 
not straightforward because of the complex PSFs which 
arise from the multiple positions and roll angles at which 
the data were taken. Here the fluxes were taken from 
the BOl catalog (converted to F = 1.2), which carefully 
calculated the total counts for the bright objects. The 
incompleteness and bias corrections were determined by 
comparing the raw number counts with those in the 6.5' 
radius region. 5 x 10^^^ erg cm^^ s~^ lies well above the 
flux at which such effects become significant. 

The 78 ks SSA22 field was analyzed using the stan- 
dard reduction tools. A complete list of sources lying in 
the S2 and S3 chips and their fiuxes were generated with 
WAVDETECT, using the assumptions of the present pa- 
per. More details can be found in M. W. Bautz et al., in 
preparation. The total area used was 0.035 square degrees. 
Only sources lying above 6 x 10"^^ erg cm~^ s~^ (above 
which incompleteness and bias corrections are small) are 
included in the present analysis. 

3. NUMBER COUNTS 

We combined the fields by forming a total list of sources 
and determining the combined incompleteness and bias- 
corrected areas sampled at each fiux. For the CDF-S we 
included the entire sample from COl, while for the SSA13 
counts from Mushotzky et al. (2000) we restricted our- 
selves to sources with fiuxes above 4 x 10~^^ erg cm~^ s^^ 
to avoid any issues of bias and incompleteness which 
might be present at lower fiuxes. The final sample con- 
sists of 373 sources with fluxes between 2 x 10~^^ and 
1.3 X 10"^"^ erg cm~^ s~^. The total corrected area ranges 
from 40 square arcminutes at the faint flux end to 0.25 
square degrees at bright fluxes. The differential counts 
are the sum of the inverse areas of the sources in the flux 



interval AS" divided by the width: 



n{S) 



E(1M.) 
AS* 



(1) 



We have normalized to a unit flux of 10^^^ erg cm~^ s^^. 
In Fig. ^ we present the differential number counts per 
square degree per unit flux (flUed squares) with la uncer- 
tainties from the Poisson error in the number of sources. 
We also show the number counts from only the CDF-N 
and Hawaii fields (open triangles) and the CDF-S counts 
from COl (open diamonds). There is a substantial dif- 
ference between these two determinations of the number 
counts at fluxes below 10~^^ erg cm~^ s~^. This dif- 
ference could reflect the different methods of analysis 
(indeed, it occurs at fluxes where modelling corrections 
are required); however, all the methods of analysis so far 
applied to the fields have shown these same differences 
(BOl; Giacconi et al. 2001b; COl), and it appears more 
probable that they result from ficld-to-field variations. 
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Fig. 2. — The combined differential counts from the four fields 
(solid boxes) with la uncertainties. The y-axis units are relative to 



a unit flux of 10 erg 



The open diamonds show the 



COl CDF-S counts, and the open triangles show the combined CDF- 
N, SSA22, and SSA13 counts. The solid lines show the power-law 
fits to the composite counts over the flux ranges above and below 
10~^^ erg cm~ s~^, while the dashed lines show the ASCA counts 
of della Cecca et al. (2000) and the BeppoSAX counts of Giommi 
et al. (2000). The dot-dash line shows Baldi et al. (2001)'s fit to 
their XMM-Newton counts. The dotted line shows COl's fit to the 
CDF-S and ASCA data. 

The number counts cannot be fit by a single power-law 
but are well fit by a broken power-law. An area weighted 
maximum likelihood fit to the counts from all four fields, 
utilizing the bias-corrected areas (Murdoch, Crawford, & 
Jauncey 1973) gives a fit above 10~^^ erg cm"^ s^^ of 



n{S) = (39 ±5) X (5/10"^'') 



-14\-2.57±0.22 



and below 10 erg cm ^ s ^ of 



n{S) = (32 ±2) X {S/IQ-''') 



-14\-1.63±0.05 



(2) 



(3) 



The normalizations are determined by matching the to- 
tal number of objects in each range, and the uncertainties 
are la. The units of n{S) are number per square degree 
per 10"^*^ erg cm^^ s^^. The two power-laws intercept at 
a fiux of 1.2 X 10^^"* erg cm~^ s~-^. A fit with the un- 
corrected areas gives a slightly steeper slope at the faint 
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end, consistent with the analytic expectation when un- 
corrected areas are used instead of bias-corrected areas 
(e.g., Schmitt and Maccacaro 1986). Fitting to only the 
CDF-N-|-Hawaii field counts gives a similar faint end slope 
(—1.61 ±0.06) but higher normalization (45 ±4), which, 
on integration, closely matches the cumulative 2 — 8 keV 
counts of BOl. COl, using only the CDF-S data and force 
fitting the bright end of the flux counts to the ASCA slope, 
inferred a faint flux slope of — 1.58±0.03 and a break point 
of 1.7 X 10~^* erg cm~^s~-^ (the dotted line in Fig. |2|). A 
systematic effect which could change the slope is the pro- 
gressive hardening of the sources as the flux drops. The 
median F in BOl drops from 1.4 at 10"^"* erg cm^^ s~^ to 
0.9 at the faintest fluxes. Using these F s to compute the 
flux conversion increases the faint end slope to —1.68. 

Cross calibration is a concern in comparing with counts 
from other satellites. Tests using a point source common 
to ASCA and Chandra give agreement at about the 10% 
level for these satellites (Barger et al. (2001)) A more de- 
tailed analysis has recently been given by Snowden et al. 
(2001) who found that the XMM and Chandra ACIS-S3 
data agree to about 5% in the 2 — 10 keV band and that 
the ASCA and Chandra fluxes agree to about 10%. The 
bright counts presented here match extremely well to the 
counts obtained by Baldi et al. (2001) using XMM-Newton 
(dash-dot line in Fig. ||) and lie only slightly below the 
ASCA counts of Ueda et al. (1999) and della Ceca et al. 
(2000) and the BeppoSAX counts of Giommi et al. (2000) 
(dashed lines in Fig. ||). 

4. DISCUSSION 

Fig. 1^ shows the contribution to the XRB versus flux. 
The peak contribution arises near the break at 1.4 x 
10^^** erg cm^^ s~^. The total contribution is 1.1 x 
10^^^ erg cm~^ s~^ deg~^ over the flux range 2 x 10^^^ 
to 10^^^ erg cm~^ s~^. Extrapolating the power-law fit 
would add a further 0.2 x 10^^^ erg cm~^ s~^ deg~^ 
above 10"^'^ erg cm~^ s^^, similar to the value ob- 
tained from the BeppoSAX number counts. Use of the 
highest ASCA counts would only increase this to 0.3 x 
10~^^ erg cm~^ s~^ deg~^. Extrapolating the faint end 
flux counts below 2 x 10"^^ erg cm^^ s^^ would add 
0.2 x 10~^^ erg cm~^ s^^ deg~^ to the total. For the 
HEAO-1 normalization the observed contribution from 
sources above 2 x 10~^^ erg cm~^ s~^ exceeds the back- 
ground, while for the BeppoSAX and ASCA normaliza- 
tions it is approximately 80% of the XRB. This would rise 
to 90% if we include the extrapolated faint end flux con- 
tribution. 



The field-to-field differences can clearly be seen 
in Fig. |, with the CDF-S and SSA22 lying below 
the CDF-N and SSA13. Between 2 x 10'^^ and 
10"^"* erg cm-2 g-i ^he CDF-N contributes 0.81 x 
10""'^"'^ erg cm~^ s~^ deg~^, while the CDF-S contributes 
0.54 X 10"-^^. This 40% difference is substantially above 
the Poisson noise expected from the number of sources. 
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Fig. 3. — The contribution to the XRB versus flux. The solid 
boxes are the measured values in the combined sample. The lines 
show the values from the power-law fits. The open boxes show the 
CDF-N, the open diamonds the CDF-S, the open upward pointing 
triangles the SSA13 field, and the open downward pointing trian- 
gles the SSA22 field. The individual fields are shown only below 
10"^* erg cm~^ s~^ where the error bars are small. 

Given these variations and the uncertainties in the flux 
calibration of both the counts and the background, we esti- 
mate the minimum contribution of the resolved sources to 
the XRB to be 1.3x10"^^ erg cm~^ s~^ deg~^. Comparing 
this number to the maximum measured value of the XRB, 
1.8x10"" erg cm"^ s^^ dcg^^ (Vecchi et al. 1999), we find 
a maximum residual of 0.5 x 10~^^ erg cm~^ s~^ deg~^. 
The true value of any diffuse component is likely to be 
considerably lower. 
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